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Abstract 



Quite good agreement has been achieved between 
computer modeling and actual performance of the 
Brookhaven 200 MeV Linac. We will present compar- 
isons between calculated and measured performance 
for beam transport through the RFQ, the 6 m trans- 
port from RFQ to the linac, and matching and trans- 
port through the linac. 

1 INTRODUCTION 

The Brookhaven 200 MeV linac serves as the injector 
for the AGS Booster and as well delivers beam to the 
Brookhaven Isotope Resource Center. It consists of a 
35 keV magnetron surface plasma source, a low energy 
beam transport (LEBT) jjj, 201 MHz radio frequency 
quadrupole (RFQ) medium energy beam transport 
(MEBT) and 200 MeV Linac §. In this last year we 
have gone through a linac upgrade to get 2.5 more 
average current (146 fx A) B. This was achieved by 
increasing repetition rate 5 to 7.5 Hz and increasing 
peak current from 25 to 39 mA. In this paper we com- 
pare computer modeling with actual performance. 




J B= 1.083 E=0.i: 




Figure 1: Ion trajectories through the LEBT and 
phase space at ion-source, middle of the line, and at 
the entrance of the RFQ. 



2 LEBT AND RFQ 

LEBT had two pulsed solenoids, two sets of x and 
y steerer, beam chopper, emittance probe, and two 
current toroids. The chopper was removed from the 
line, making the line shorter by 70 cm. Computer 
modeling of this line showed that we should move the 
1st solenoid as close to the ion source as possible to 
reduce the beam size in the 1st solenoid and the second 
solenoid as close to the RFQ as possible to increase the 
convergence angle required by the RFQ acceptance. 
Figure 1 shows the ion trajectories through this line 
and phase space at the exit of the ion source, middle of 
the line and entrance of RFQ. The RFQ acceptance is 
shown as a solid line ellipse. This calculation assumed 
that the beam space charge is neutralized. 

Shortening of the line resulted in lower measured 
emittance. Due to lower emittance and better match- 
ing transmission through the RFQ was improved by 
about 10 percent. Figure 2 shows the transmission as 
a function of input beam current. 
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Figure 2: RFQ transmission as a function of input 
current. The solid line indicates the simulation and 
crosses the measured values. 



*Work performed under the auspices of the U. S. Department 
of Energy. 
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Figure 3: MEBT, 750 keV Transport Line. 

3 MEBT AND LINAC 

This is where we have recovered most of the beam 
losses. This line is 6 meter long and is shown in Fig- 
ure 3. It consists of four triplets, three bunchers, one 
slow chopper, one fast chopper, two emittance mea- 
surement units, three current transformers, two sets 
of x, y steerers, and a dipole to accomodate polarized 
beam coming at 60 deg angle. The ideal match be- 
tween RFQ and DTL could have been obtained with 
a 5 j3\ long FODO lattice with quadrupole spacing 
about [3X and at least two bunchers. But the require- 
ment of beam chopping and polarized beam dictated 
a triplet solution ||. 

The first quadrupole after the RFQ was too far; by 
the time the beam reached the quadrupole it had gone 
through a waist in the x plane, hence was diverging 
in both planes. No matter which polarity quad one 
puts, beam size in the other direction is very big. Also 
the longitudinal beam size is too big before it reaches 
the first buncher. To improve the capture and trans- 
mission of the beam in MEBT, the RFQ end flange 
at the high energy end was modified to accommodate 
a permanent magnet quadrupole (PMQ). The PMQ 
was similar to one as used in the SSC DTL §. We 
have also rearranged the gate valve and current trans- 
formers at the beginning of the line, and also measured 
and aligned all the quadrupoles very carefully. Mea- 
surement as well as simulation showed that as little 
as a 1.5 degree quadrupole rotation can increase the 
emittance by 50 %. The last quadruplet was changed 
to a triplet to reduce coupling. Figure 4, shows the 
measured and calculated phase spaces after the sec- 
ond buncher. Table 1 shows calculated (TRACE3D) 
and measured Twiss parameters for Figure 4. 

Table 2, compares measurements and PARMILA re- 
sults at various locations. We believe that lower values 
for currents after the second buncher and and at the 
entrance of the Linac are caused by the grids in the 



Figure 4: Measured (upper two) and calculated (lower) 
phase space after the second buncher. 



Planes 
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0v 




Meas. 


-0.60 


0.79 


9.71 


-1.59 


1.61 


13.82 


Cal. 


-0.31 


0.66 


9.80 


-1.79 


1.67 


14.00 



Table 1: Calculated and measured Twiss parameters 
after the second buncher. (3 is in mm/mrad and e 
(minor., RMS) in 7r mm mrad. 

buncher drift tubes (four in each buncher). These grids 
are placed in the bunchers to reduce RF defocusing ef- 
fects. Emittance measurements at 200 MeV are done 
using profiles at five places. Agreement between cal- 
culations and measured Twiss parameters at 200 MeV 
is poor because there are 295 quadrupoles, and cali- 
bration and misalignment errors are not known to a 
good accuracy. 



Location 


94-95 


95-96 


Sim. 


Meas. 


Sim. 


Meas. 


Current mA 


RFQ 


50.0 


50.0 


62.9 


62.9 


Buncher 


50.0 


41.0 


62.9 


57.8 


MEBT 


42.4 


37.3 


62.9 


53.2 


Tank 1 


27.4 


28.4 


37.1 


37.7 


Tank 9 


26.2 


26.7 


36.4 


35.9 


Emittance, (nor, RMS) 7r mm mrad 


RFQ 


0.400 


0.400 


0.375 


0.375 


Buncher 


0.44 


0.56 


0.47 


0.57 


200 MeV 


1.32 


2.8 


1.85 


1.92 



Table 2: Comparison between simulations and mea- 
sured beam parameters. 



4 ALGEBRAIC 
RECONSTRUCTION 
TECHNIQUE (ART) 

A radiograph of the beam at the BLIP target taken last 
year showed a tilted ellipse in the x-y plane. Sources of 
this coupling can be quad rotation or vertical offset in 
the dipole. This triggered the need for an x-y density 
profile. We found that algebraic reconstruction tech- 
nique (ART) could help us. ART was introduced by 
Gordan, Bender and Herman Q for solving the prob- 
lem of three dimensional reconstruction from projec- 
tions. The ART algorithms have a simple intuitive ba- 
sis. Each projected density is thrown back across the 
reconstruction space in which the densities are itera- 
tively modified to bring each reconstructed projection 
into agreement with the measured projection. The re- 
construction space is an n x n array of small pixels, p 
is grayness or density number which is uniform within 
the pixel but different from other pixels. Assume P 
is a matrix of m x n 2 and the m component column 
vector R. Let pi.j denote the (i,j)th element of P , and 
Ri demote the ith element of reconstructed projection 
vector R. For 1 < i < m, Nj is number of pixels under 
projection R^, defined as N t — J2jPl.j- The density 
number denotes the value of pj after q iterations. 

After q iterations the intensity of the ith recon- 
structed projection ray is 



and the density in each pixel is 
Ri - R? 



Pj = pj+pi,j- 



with starting value = 



where R^ is the measured projection and, 
m,if (q+1) is divisible m 

the remainder of dividing (q+l)by m, otherwise 



and, 



0. 



p q ={ PT 



if < 
if < p~ q < 1 

i; if pj q > i 

It is necessary to determine when an iterative algo- 
rithm has converged to a solution which is optimal ac- 
cording to some criterion. We are using as the criteria 
of convergence the discrepancy between the measured 
and calculated projection elements 



D" = 



1 ^(Ri-R 



q\ 2 



N, 



We have added a third wire at 45 degree in wire 
scanners in two places. Figure 5 compares the mea- 
sured and reconstructed profiles in the BLIP transfer 
line |8| after the 1st octupole and figure 6 shows the 
reconstructed 3D density distribution. 
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Figure 5: Beam projection on x, y, and 45 degree 
planes. 
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Figure 6: Reconstructed 3D density distribution us- 
ing ART. 
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